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Abstract Distance measurements of gamma-ray pulsars are challenging questions in present 
pulsar studies. The Large Area Telescope (LAT) aboard the Fermi gamma-ray observatory 
discovered more than 70 gamma-ray pulsars including 24 new gamma-selected pulsars which 
nearly have no distance information. We study the relation between gamma-ray emission effi- 
ciency (r) — Lry/E) and pulsar parameters for young radio-selected gamma-ray pulsars with 
known distance information in the first gamma-ray pulsar catalog reported by Fermi/LAT. We 
have introduced three generation order parameters to describe gamma-ray emission properties 
of pulsars, and find the strong correlation of r\ — £3 a generation order parameter which reflects 
7-ray photon generations in pair cascade processes induced by magnetic field absorption in 
pulsar magnetosphere. A good correlation of r\ — £>lc the magnetic field at the light cylinder 
radius is also found. These correlations would be the distance indicators in gamma-ray pulsars 
to evaluate distances for gamma-selected pulsars. Distances of 25 gamma-selected pulsars are 
estimated, which could be tested by other distance measurement methods. Physical origin of 
the correlations may be also interesting for pulsar studies. 
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1 INTRODUCTION 

Before 2008, only 7 gamma-ray pulsars are known in nature (Thompson 2001). The launch of the Fermi 
Gamma-ray Space Observatory in June 2008 completely changed the status in studies of gamma-ray pul- 
sars. The first published catalog of gamma-ray pulsars (Abdo et al. 2010) contains 46 gamma-ray pulsars 
including 8 millisecond pulsars, 21 young radio pulsars and 17 gamma-selected pulsars. After more than 
one and half years of all-sky survey observations by Fermi/LAT, more than 70 gamma-ray pulsars were 
discovered, including 25 gamma-selected pulsars (see reviews by Ray & Saz Parkinson 2010). High sensi- 
tivity of the Fermi/LAT makes a new era for pulsar discoveries, specially for the population of radio-quiet 
gamma-ray pulsars. 

The distance measurement of pulsars is always a difficult problem in pulsar studies. Trigonometric 
parallax measurements of radio pulsars are the reliable method, but are only available for the nearby pulsars 
(< 0.4 kpc) specially for a few radio millisecond pulsars (e.g. Lommen et al. 2006). The most common way 
to obtain radio pulsar distance is based on the computation from dispersion measurement (DM) coupled 
to an electron density distribution model like NE 2001 model (Cordes & Lazio 2002), which have been 
applied to most radio pulsars (e.g., Johnston et al. 1996; Keith et al. 2008). The pulsar distance can be also 
estimated from kinematic model: the distance of possible associated objects (supernova remnants, pulsar 
wind nebulae, star clusters, or HII regions) could be measured from Doppler shift of absorption or emission 
lines in HI spectrum together with the rotation curve model of the Galaxy (e.g., Robert et al. 1993; Camilo 
et al. 2006). The distance of some pulsars with X-ray emissions can be estimated from X-ray observations 
of the absorbing column (e.g. Romoni et al. 2005) or from correlations in X-ray luminosities versus spin- 
down power or photon index (Becker & Tmemper 1997; Possenti et al. 2002; Gotthelf 2003; Wang 2009; 
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and references therein). These methods may be available for radio or even X-ray pulsars, but for gamma- 
selected pulsars if no possible associated objects, we would have no any information on their distance. 

It is well known that X-ray luminosity has the correlation with pulsar's spin-down power: L x oc E in 
soft X-ray bands (0.1 - 2.4 keV, Becker & Truemper 1997), and L x oc E z / 2 in hard X-ray bands (> 2 
keV, see Saito 1998; Cheng et al. 2004; Wang 2009). Based on the EGRET pulsars, Thompson et al. (1999) 
found a possible correlation of L 7 oc E 1 / 2 . For the larger sample of gamma-ray pulsars in Abdo et al. 
(2010), the young pulsars looks to still follow this relation with a large scattering factors of more than 10 
but millisecond pulsars follow a different relation (see Fig. 6 of Abdo et al. 2010). This correlation was used 
to estimate some gamma-selected pulsars (Saz Parkinson et al. 2010). Moreover, the relation of L 1 oc E 1 / 2 
may not be intrinsic, for the young gamma-ray pulsars in Fig. 6 of Abdo et al. (2010), we find a fitting 
function of L 7 oc E 0,7 . 

Gamma-ray emission efficiency (rj = L 7 /E) is an important parameter in gamma-ray pulsars, which 
varies for different populations of pulsars. In this paper we study the relations of gamma-ray emission effi- 
ciency versus some pulsar parameters: spin period, age, magnetic field at light cylinder, and three generation 
order parameters. We will show results of these relations and good correlations would be pulsar distance 
indicators for gamma-selected pulsars. 

2 GAMMA-RAY EMISSION EFFICIENCY VERSUS PULSAR PARAMETERS 

Gamma-ray emission efficiency is defined as -q — L 7 /E, where the spin-down power E = 4ir 2 IPP~ 3 
taking / = 10 45 g cm 2 , P is the period of pulsar in units of second. L 1 = And 2 /qF 7 , where F 1 is the 
gamma-ray flux above 100 MeV detected by Fermi/LAT. The radiation open angle factor fa is model- 
dependent, and may depend on the magnetic inclination and observer angles, which could be obtained 
using pulse profile information (e.g., fn^l for 8 gamma-ray pulsars estimated by Watters et al. 2009). For 
simplicity, we use /o = 1 in this paper similar to Abdo et al. (2010). In the gamma-ray pulsar catalog given 
by Abdo et al. (2010), 21 radio-selected young pulsars and 7 gamma-selected pulsars have the distance 
measurements or estimation. However, gamma-ray emission efficiency of some gamma-selected pulsars 
are higher than 1, the maximum radiation efficiency in physics, suggesting overestimation of the distance 
for some gamma-selected pulsars. Millisecond pulsars may have different properties from young pulsars, so 
we do not consider 8 millisecond pulsars in the catalog here. Finally, we use these 21 radio-selected pulsars 
for the analysis in this section. The efficiency rj distributes from 0.1% (like Crab pulsar) to near 100%. 

We will first show the relations between rj versus three well known pulsar parameters: P, r and £>lc- 
Most importantly, we have introduced the generation order parameters for pulsars (see details in Wang & 
Zhao 2004 and therein references) which can be used to describe gamma-ray properties in pulsars. In §2.2, 
the relations between i] versus three generation order parameters will be studied. 

2.1 r) versus P, t and Blc 

In Fig. 1, we plot diagrams of rj versus P, r and B^c for 21 young gamma-ray pulsars, respectively. 
t = P/P is the pulsar's characteristic age, £?lc ~ 2.94 x 10 8 (PP~ 5 ) 1 / 2 is the magnetic field at the light 
cylinder (R LC = cP/2tt). 

In the diagram of 77 — P, the data points of spin period are scattering, no significant correlation is found. 
But rj seems to have the correlations with the other two pulsar parameters: age and the magnetic field at the 
light cylinder. The linear function is used to fit the correlations (solid lines in Fig. 1): 

logr? = -(4.73 ± 0.31) + (0.68 ± 0.08) logr (1) 
with a standard deviation of a ~ 1.63 and the probability value (p-value for t-test) of 1.09 x 10~ 4 ; 

log rj = (2.23 ± 0.32) - (0.88 ± 0.10) logS LC (2) 

with a <~ 1.56, and a p-value of 2.90 x 10 -5 . 

The gamma-ray efficiency generally become higher with older age and smaller £?lc- From the evalua- 
tion of the standard deviation values a and p-values, the relation of 77 — £?lc is the better one. 




Fig. 1 Gamma-ray emission efficiency 77 of 21 young gamma-ray pulsars versus three pulsar 
parameters: spin period P, age r and the magnetic field at light cylinder Plc- V shows the 
correlations with two pulsar parameters t and -Blo an d the solid lines display the best fitting 
functions. See the text for details. 




Fig. 2 Gamma-ray emission efficiency 77 of 21 young pulsars versus three generation order pa- 
rameters (i_3 in gamma-ray pulsars, rj has no significant correlation with £1 but has the corre- 
lation with (2 and £3, suggesting that the magnetic field dominates the gamma-ray absorption in 
cascade processes. The sold lines show the best fitting function. See the text for details. 



2.2 7] versus generation order parameters 

The concept of generation is provided to describe pair cascade processes in gamma-ray pulsars (Zhao 
et al. 1989, Lu & Shi 1990). Based on Ruderman-Sutherland scenario (Ruderman & Sutherland 1975), 
passing through the polar gap, e + /e~ are accelerated to a high energy with typical Lorentz factor = 

6.0 x 10 7 P 1 / 14 P 1 ^ 1 ^ 14 , where P15 the derivative in units of 10~ 15 s s -1 . These first generation particles 
will move along the curved magnetic field lines and emit high energy curvature radiation (the first generation 
photons) with photon energy typically at 

Ei = ~ll « 3.2 x 10 10 F- 2 / 7 Pr 5 3/14 eV, (3) 
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where R c w 1.8 x 10 7 P 1 / 2 cm is the curvature radius of field line here. These primary photons could be 
converted into secondary e+/e~ pairs in both open and closed magnetic field line regions near the neutron 
star surface due to magnetic pair creation (Halpern & Ruderman 1993). And the condition of these photons 
to create e + /e~ pairs is (Sturrock 1971; Ruderman & Sutherland 1975) 

S i ~ 1 (4) 



2m e c 2 B c 15 ' 

where B(r s ) is the local magnetic field at the position of r s , and B c = m^c 3 jeh = 4.14 x 10 13 G is the 
critical magnetic field. These e + /e~ can emit the second generation photons through synchrotron radiation 
with a characteristic energy E 2 . If E 2 is high enough, the further e + /e~ pairs (the third generation) can be 
produced under the condition similar to Eq. (4), 

E 2 B(r s ) _ xo 5 



2m e c 2 B c 15 ' 

where Xo/15 ~ 1/9 — 1/12 (Sturrock 1971). Then pair cascade processes occur. 

Concerning this idea, Lu et al. (1994) introduced the generation order parameter (GOP) to characterize 
a pulsar. They considered the conversion of high energy photons into e + /e~ pairs through electric fields, 
and defined the first GOP as 

c =1 , l-(ll/7)logP + (4/7)logA 5 (6) 
1 3.56 - logP - logAs 

Wei et al. (1997) considered absorption of high energy photons by the effect of both magnetic and 
electric fields, define the second GOP as, 

, i ^ 0.8-(2/7)logP+(2/7)logA 5 n , 
C 2 = l + U ■ (7) 

The concept of generation was initially considered in the scheme that the 7-ray photons is absorbed and 
conversed into e + /e~ through only magnetic fields (Zhao et al. 1989), so we defined the third GOP based 
on the magnetic field absorption effects as (Wang & Zhao 2004) 

, 1 , 0.6-(ll/14)logP+(2/7)logA 5 m 
Ca = l + u ■ (8) 

GOPs are used to describe cascade processes and characterize the spectral properties of pulsars. If a 
pulsar can emit gamma-rays, the GOPs must be larger than 1 (i.e., the first generation gamma-ray photons 
must exist). In addition, the GOPs are proved to be correlated with the gamma-ray photon index: softer 
gamma-ray photons with larger GOPs based on the EGRET pulsar sample (Lu et al. 1994; Wei et al. 1997). 
Then according to the definition of GOPs, the first generation pairs emit high energy gamma-rays (i.e., 
> 100 MeV), with larger GOPs, more first generation pairs are transferred into next generation pairs with 
lower energy which emit more soft gamma-rays and X-rays. So given a total emission rate, efficiency to 
GeV gamma-rays (77) becomes lower with larger GOPs. 

In Fig. 2, we plot the diagrams of 77 versus three GOPs (C1-3) respectively. 77 has no correlation with 
Ci but has correlation with other two GOPs ( 2 , (3, implying that magnetic fields dominate the absorption in 
pair cascade processes, consistent with our previous results ( Wang & Zhao 2004). These correlations also 
suggest that GOPs (( 2 , £3) can describe gamma-ray properties of pulsars. In Fig. 2, the solid lines show the 
best fitting functions for the relations of 77 — ( 2 and 77 — £3: 

log r\ = (4.98 ± 0.45) - (3.00 ± 0.21)C 2 (9) 
with a ~ 1.69 and ap-value of 1.94 x 10~ 4 ; 

log r\ = (5.49 ± 0.24) - (2.86 ± 0.11)Cs (10) 

with a ~ 1.34 and a p-value of 1.01 x 10~ 6 . The correlation between 77 — ( 3 is stronger with the smaller 
standard deviation and p- value. 
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Table 1 The estimated distances of 25 gamma-selected pulsars. di_2 denotes the distance range 
calculated from relations of 77 — C3 and 77 — £?lc> respectively. d 3 is the estimated distance from 
other methods with references provided. 
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3 POSSIBLE DISTANCE INDICATORS FOR GAMMA-RAY SELECTED PULSARS 

In §2, the relations between r\ and six pulsar parameters: P, r, i?LC and three GOPs £1, ^2, C3 are stud- 
ied. From the values of standard deviation and p-values after fittings, the correlations of r\ — £3 is stronger 
than others, and the correlation of rj — B^c could also be acceptable. In this paper we do not consider the 
physical origin in these correlations. These pulsar parameters can be estimated by two fundamental mea- 
surement parameters P and P which are relatively easily observed. The gamma-ray emission efficiency is 
sensitively dependent on distance measurement which is very difficult at present, specially nearly impossi- 
ble for gamma-selected pulsars. With the strong correlation of 77 — (3, we have a possible way to estimate a 
reliable distance for gamma-ray pulsars with only known P, P and F 7 . 

In the catalog of Abdo et al. (2010), 17 gamma-selected pulsars are listed and most of them have no 
any distance information. Saz Parkinson et al. (2010) claimed detections of 8 new gamma-selected pulsars 
in blind frequency searches of Fermi LAT data. In Table 1, we use the distance indicator obtained by 77 — (3 
correlation to estimate the distances of the 25 gamma-selected pulsars. For comparison, we also give the 
predicted distance values calculated using the relation of rj — B^q. From Table 1, we find the evaluated 
distances {d\, cfe) by r\ — £3 and 77 — £>lc correlations are similar, suggesting that these two distance 
indicators can be checked by each other. 

In Table 1, we also collected the distance information (d^) for some gamma-selected pulsars from other 
measurements or estimations. For the Geminga pulsar, we estimate the distance of 0.19 ± 0.07 kpc which 
is well consistent with the distance value of 0.25l° kpc from the optical trigonometric parallax measure- 
ment (Faherty et al. 2007). For PSR J1836+5925, we estimate its distance as <~ 0.3 kpc (corresponding 
to an efficiency of ~ 55%) which is also well below the upper limits of 0.8 kpc according to its thermal 
X-ray spectrum (Halpern et al. 2007). For other gamma-selected pulsars, our estimated distance values are 
generally below those from other methods, but may be more reliable. According to the distance estimated 
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Fig. 3 Gamma-ray emission efficiency 77 versus £>lc an d C3 f° r both 21 young gamma-ray pul- 
sars (solid circles) and 8 millisecond gamma-ray pulsars (open circles). Millisecond pulsars still 
generally follow the relations of 77 — _Blc an d V — C3 m young pulsars, but they may have a nearly 
constant gamma-ray radiation efficiency of 77 ~ 10%. 



from the 77 — C3 relation, the gamma-ray efficiency 77 is general below 1 . The estimated efficiency of PSR 
J2021+4026 is about ~ 0.16 (corresponding to d ~ 0.4 kpc) according to the 77 — ^3 relation, compared 
to 77 ~ 0.9 — 3.6 (corresponding to distance of 1 - 2 kpc) from kinematic model method on the possible 
association (Landeckeret al. 1980). 



4 SUMMARY AND DISCUSSION 

In this article, we studied the possible correlations between gamma-ray emission efficiency 77 and 6 pulsar 
parameters: P, r, Blc an d three generation order parameters (is using 21 young radio-selected gamma- 
ray pulsars in Abdo et al. (2010). We find the strong correlation between 77 — (3. Based on the concept of 
the GOPs, lager GOPs imply that more high energy photons are transferred to softer photons (X-rays). The 
good correlation of 7/— £3 suggests that GOPs can describe gamma-ray emission properties of young pulsars, 
and the magnetic field absorption effects dominate in pair cascade processes in pulsar magnetosphere. 
This intrinsic correlation can be used to estimate distances for gamma-selected pulsars which have no any 
distance information yet. The good correlation of 77 — £>lc is a l so found, which can be also used as the 
other distance indicator in gamma-ray pulsars for a double check. 

Millisecond pulsars (MSPs) have not be included in our analysis, though their distances are generally 
measured by optical trigonometric parallax and DM methods. MSPs with much smaller P and P have a 
much older characteristic age (r ~ 10 9 yr). The values of (1 and (2 are below 1 or near 1, making MSPs 
as non-gamma pulsars if these two GOPS are still applicable to MSPs. However, in parameter spaces of 
£?lc an d £3, MSPs are similar to young pulsars. In Fig. 3, we plot the diagrams of 77 — Blc an d 77 — ^3 
including both 21 young pulsars and 8 MSPs in the first gamma-ray pulsar catalog (Abdo et al. 2010). 
MSPs seems to still follow the behaviors of young pulsars: higher efficiency with smaller values of B^c 
and £3. In the same time, gamma-ray emission efficiency of MSPs could also be thought to keep constant 
77 ~ 10% (also see Fig. 6 of Abdo et al. 2010). So MSPs may have different gamma-ray emission properties 
from young pulsars, like multi-pole magnetic field assumption in MSPs (Ruderman 1991; Zhang & Cheng 
2003), or different emission open angles (taken as /n ~ 0.5, Fierro et al. 1995). Present discoveries of 
MSPs are generally done through radio timing, and the blind search for MSPs by Fermi/LAT is a very 
important project in future, but quite difficult at present specially for MSPs in binaries. Then the distance 
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Fig. 4 The distance distributions of three classes of gamma-ray pulsars: gamma-selected pul- 
sars, radio-selected pulsars and millisecond pulsars. The distances of gamma-selected pulsars 
are taken from the column d\ of Table 1 according to the distance indicator of the r\ — (3 corre- 
lation. 



indicators of 77 — £?lc an d V — C3 relations could be the secondary way for distance information of MSPs 
after trigonometric parallax or DM methods. 

The GOP model was originally proposed based on the polar-cap accelerator scenario. The present 
Fermi/LAT may support that gamma-ray emission in pulsars mainly comes from the spatially extended 
regions reaching a good fraction of the light-cylinder radius (e.g., Abdo et al. 2010). The production of 
the secondary pairs in polar-cap activity is also different from that in the scenarios of outer-gap models 
or slot-gap models (e.g., Cheng, K. S. et al. 2000; Muslimov & Harding 2004). Then the new model of 
generation order parameters may be developed in the extended regions from the polar-cap regions to near 
the light-cylinder radius. This GOP model would be more complicated but could be considered in the next 
work. Anyway, the correlations of 77 — £3, 77 — £?lc for gamma-ray pulsars suggest that the gamma-ray 
luminosity may depend on two fundamental pulsar parameters P and P. The function of P and P could 
well predict gamma-ray emission luminosity, which can be used to trace the distance of gamma-ray pulsars. 

Fig. 4 shows the distance distributions of three classes of gamma-ray pulsars: gamma-selected pulsars, 
radio-selected pulsars and millisecond pulsars. The distances of gamma-selected pulsars are provided by 
the distance indicator of the 77 — £3 relation (see Table 1). Gamma-ray loud millisecond pulsars distribute 
at a distance peak around 0.3 kpc because MSPs generally have lower spin-down powers. Gamma-selected 
young pulsars distribute at the distance peak of ~ 1.2 kpc, while radio-selected young pulsars distribute 
at the distance peak of ~ 2.5 kpc. This difference in distance distributions for two classes of gamma-ray 
young pulsars may involve further interest. The nearby unresolved radio-quiet gamma-ray pulsars could 
contribute to diffuse gamma-ray background specially for the high-latitude pulsars located in the Gould 
Belt (Wang et al. 2005). 

Before the Fermi era, only one gamma-selected pulsar Geminga was known. Now 25 gamma-selected 
pulsars are discovered, greatly improving our knowledge of gamma-ray pulsar family. Much more gamma- 
selected pulsars would be detected by future deeper sky surveys of Fermi/LAT. The distance indicators 
presented in this paper will give the distance information for gamma-selected pulsars, which will be helpful 
for study in gamma-ray emission properties of this pulsar population. It is still hopefully expected that more 
gamma-ray pulsars (young) have trigonometric parallax measurements or more precise DM model, which 
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can check the validity of the distance indicators (i.e., 77 — ^3, 77 — £?lcX an d improve the distance indicators 
in gamma-ray pulsars. 
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